Collective excitations on a surface of topological insulator 
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We study collective excitations in a helical electron liquid on a surface of three-dimensional topo- 
logical insulator. Electron in helical liquid obeys Dirac-like equation for massless particless and 
direction of its spin is strictly determined by its momentum. Due to this spin-momentum locking, 
collective excitations in the system manifest themselves as coupled charge- and spin-density waves. 
We develop quantum field-theoretical description of spin-plasmons in helical liquid and study their 
properties and internal structure. Value of spin polarization arising in the system with excited 
spin-plasmons is calculated. We also consider the scattering of spin-plasmons on magnetic and non- 
magnetic impurities and external potentials, and show that the scattering occurs mainly into two 
side lobes. Analogies with Dirac electron gas in graphene are discussed. 
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I. INTRODUCTION 

Topological insulator is a new class of solids with non- 
trivial topology, intrinsic to its band structure. Theo- 
retical and experimental studies of topological insulators 
grow very rapidly in recent years (see [1, 2] and references 
therein). Three-dimensional topological insulators are in- 
sulating in the bulk, but have gapless surface states with 
numerous unusual properties. These states are topolog- 
ical^ protected against time-reversal invariant disorder. 
When gap is opened in surface states by time-reversal or 
gauge symmetry breaking, a spectacular magnetoelectric 
effect arises [3, 4]. 

Recently a "new generation" of three-dimensional 
topological insulators (the binary compounds Bi2Se3, 
Bi2Tc3 and other materials), retaining topologically pro- 
tected behavior at room temperatures, were predicted 
and studied experimentally [5-7]. Band structure of the 
surface states of these materials contain a single Dirac 
cone, where electrons obey two-dimensional Dirac equa- 
tion for massless particles. Direction of electron spin in 
these states is strictly determined by their momentum, 
so these states can be called as "helical" ones. Surface 
of topological insulator can be chemically doped, form- 
ing charged helical liquid. The spin-momentum lock- 
ing leads to interesting transport phenomena including 
coupled diffusion of spin and charge [8], inverse galvano- 
magnetic effect (generation of spin polarization by elec- 
tric current) [9] and giant spin rotation on an interface 
between normal metal and topological insulator [10]. The 
spin-momentum locking offers numerical opportunities 
for various spintronic applications. 

Collective excitations (plasmons) in helical liquid on 
the surface of topological insulator was considered in 
[11]. It was shown that due to spin-momentum lock- 
ing responses of charge and spin densities to external 
electromagnetic field are coupled to each other. Thcrc- 
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fore the plasmons in the system should manifest them- 
selves as coupled charge- and spin-density waves and can 
be called "spin-plasmons". In [12] application of spin- 
plasmons in spin accumulator device was proposed. Also 
surface plasmon-polaritons under conditions of topologi- 
cal magnetoelectric effect were considered in [13]. 

Properties of the states on a surface of three- 
dimensional topological insulator are similar to those 
of electrons in graphene. Graphene is unique two- 
dimensional carbon material with extraordinary elec- 
tronic properties [14-16]. Its band structure contains two 
Dirac cones with electrons behaving as massless Dirac 
particles in their vicinities. Graphene is a perspective 
material for nanoelectronics due to large carrier mobil- 
ities at room temperature. Electronic interactions and 
collective excitations in graphene have been extensively 
studied (see [17] and references therein). In particular, 
the properties of plasmons [18, 19] and hybrid plasmon- 
photon [20] and plasmon-phonon [21, 22] modes were in- 
vestigated theoretically and experimentally. It was real- 
ized recently that graphene is a fertile ground for quan- 
tum plasmonics [23] due to very small plasmon damping. 

In this article, we develop quantum field-theoretical 
formalism to describe plasmons in graphene and spin- 
plasmons of on a surface of three-dimensional topological 
insulator based on random phase approximation. Prob- 
lems of excitation, manipulation, scattering and detec- 
tion of single plasmons can be conveniently considered 
using this approach. Thus, this approach can be espe- 
cially useful for problems of plasmon quantum optics and 
quantum plasmonics. We use our approach to study in- 
ternal structure and properties of spin-plasmons in a he- 
lical liquid. 

The rest of this paper is organized as follows. In section 
2, we present a brief description of electronic states on a 
surface of topological insulator and in graphene. In sec- 
tion 3, quantum field-theoretical description of plasmons 
in Dirac electron gas is given. In section 4, we consider 
internal structure of plasmons. In section 5, different 
consequences of spin-momentum locking are considered. 
Scattering of plasmons on impurities and external poten- 
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tials is considered in section 6, and section 7 is devoted 
to conclusions. 



II. DIRAC ELECTRONS 

The low-energy effective Hamiltonians for electrons in 
helical liquid [24] and graphene [16] are 

helical liquid :H = v F (p x a y - p y a x ), (1) 
graphene:i/ = v F (p x a x + p y a y ), (2) 

where the Fermi velocities of electrons v F are 6.2 x 10 m/s 
for the topological insulator Bi2Se3 and 10 6 m/s for 
graphene; the Pauli matrices a x and a v act in the spaces 
of electron spin projections (helical liquid) or sublattices 
(graphene). The eigenfunctions of (l)-(2) can be written 
as e lp r |/ P7 )/v / S', where S is the system area and |/ P7 ) 
is the spinor part of the eigcnfunction, corresponding to 
electron with momentum p from conduction (7 = +1) or 
valence (7 = —1) band: 



helical liquid: |/ p7 ) 
graphene: I / P7 ) 





' e -iVp/2 















(3) 
(4) 



(here and below we assume K = 1). Another difference 
between helical liquid on a surface of topological insula- 
tor and electron liquid in graphene is additional fourfold 
degeneracy g = 4 of electronic states in graphene by two 
spin projections and two nonequivalent valleys. 

The value of electron spin in helical liquid and of pseu- 
dospin graphene in the state |/ P7 ) is 

helical liquid: (/ P7 |<r|/ P7 ) =7[zxp], (5) 
graphene: (/ P7 1 er | / p7 ) = 7 p, (6) 

where p and z are unit vectors directed along the momen- 
tum p and the z-axis respectively. We see that, in helical 
liquid, the spin of electron lies in the system plane and 
makes an angle 90° (in counterclockwise direction in the 
conduction band and inversely in the valence band) with 
its momentum. In graphene, the sublattice pseudospin 
of electron is directed along its momentum in conduction 
band and opposite to it in the valence band. Physically, 
a definite direction of the pseudospin in the system plane 
corresponds to definite phase shift between electron wave 
functions on different sublattices. 

A starting point for quantum field-theoretical consider- 
ation of plasmons on the surface of topological insulator 
and in graphene is the many-body Hamiltonian of elec- 
trons with Coulomb interaction between them: 

h = E ^7 a p7 a P7 + 25 E E E v i 

P7 qpp' 7172 7(7^ 



X (/p+q,7j l/p7i) \/p'— q,72 l/p'72) 
Xa p+q,7; a P '-q,7^ a P'72 a P7l! 



(7) 



where the destruction operator for electron with 

momentum p from the band 7, £ P7 = ^v-pp — p is its ki- 
netic energy measured from the chemical potential p, and 
V q = 2ire 2 /eq is the two-dimensional fourier transform of 
Coulomb interaction potential screened by surrounding 
three-dimensional medium with a dielectric susceptibil- 
ity £. 



III. DESCRIPTION OF PLASMONS 

To investigate the properties of plasmons in Dirac elec- 
tron gas, we use equation-of-motion approach [25]. We 
treat a plasmon as a composite quasiparticle, consisting 
of electron-hole pairs with the common total momentum 
q. Thus, the creation operator for plasmon with momen- 
tum q can be written as: 



)+ 



E 

P77' 



Pk U p+q,Y U PT 



(8) 



Here the coefficients C^ 7 are the weights of intraband 
(7 = 7') and interband (7 7^ 7') single-particle transi- 
tions, contributing to the wave function of plasmon. 

The plasmon creation operator should obey Hciscnberg 
equation of motion 



[H,Q+] = n q Q+ 



(9) 



where f2 g is plasmon frequency. We start from equation 
of motion for single electron-hole pair, which can be de- 
rived using (7): 



(e 



p+q>7' 



£p7) a P 4 



5P7^p+q,7' a P7 



Hi a p+q,7' a P7 



q' 



( ^P+q-q',7" l/p+q,7') a P +q-q' j7 " a P7 



7" 



" (/p7 I /p+q' ,7" ) a p+q,7' a p+q' ,7' 



(10) 



where we have introduced the Fourier transform of the 
density operator: 



Pt = E^p+i.7'I/i 



P7/°p+q,7' a P7- 



(11) 



The right-hand side of the equation (10) contains prod- 
ucts of four fermionic operators. To reduce it to prod- 
ucts of two fermionic operators we use the random phase 
approximation (RPA); its applicability will be discussed 
below. According to RPA [25], the operator products in 
the last two lines of (10) can be replaced by their average 
values in ground state |0) of the system: 



a p'7 ,a P7 



l a P '7' a P7 



|0) = 6 PP .6. 



pp' U*yry' fipy 5 



(12) 



where the occupation number for electronic states 

with momentum p from the band 7. For electron-doped 
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Dirac liquid at T = (see also the remark at the end 
of this section), we have n p+ = Q(pp — |p|), n p _ = 1, 
where p-p = ji/vp is the Fermi momentum. In the case of 
hole doping, all characteristics of plasmons are the same 
due to electron-hole symmetry. 

Thus, the equation of motion (10) for electron-hole pair 
with taking into account the RPA assumption (12) takes 
the form: 



Hi a p+q, 7 ' a P7 
V„ 



£ P7 )a 



P7/ ,u p+q,7 /U 'P7 



g Pq ifpilfp+m') ( n P7 



(13) 



Combining the definition of plasmon creation operator 
(8) with the equation of motion for plasmon (9) and single 
electron- hole pair (13), we obtain the system of equations 



for the coefficients C^ 7 : 



Va 



IPI/ 



(fi, + Cp7 Cp+q,7')Cpq 7 — "c"(/p+q>7'l/) 

X E Cp'q(/p'r|/p'+q.r') (lip'r - V+q-r') ■ (14) 



Introducing infinitesimal damping into the plasmon 
frequency and denoting 

ff, = ^ Cpq </pr|/p+q,r') (n pT - n p+q , T ,) , (15) 



prT' 



we can find the plasmon wave function from (14) as 
(Jp+q,7'l/p7)-^q 



rn 7 
°pq 



£p7 ?p+q,' 



id' 



(16) 



The normalization factor iV q can be determined from 
the commutation relation for plasmon operators, which 
should be satisfied on the average in the ground state: 



<0|[Q q ,Q+]|0) = «W- 
Substituting (8) in (17), we get 

E D 7'7 = 1> 

77' 



(17) 



(18) 



where the total weights of intraband (D ++ ) and inter- 
band (D_| + D |_ = 1 — D ++ ) electron transitions, con- 
tributing to the plasmon wave function (16), are: 



^'7=E| C Pq 7 ( 



(19) 



To find the plasmon frequency Q q , we can substitute 
(16) in (15), arriving to the standard RPA equation (see 
it for the helical liquid [11] and for graphene [18, 19]): 



v q u( q ,n q 



0. 



(20) 



where the polarization operator of Dirac electron gas is 
introduced: 



n(<?,w) = K/p+q,7'l/p7> 

P77' 



W + £p7 - Cp+q-7' + iS 



(21) 



The degeneracy factor g is 1 for helical liquid on the 
surface of topological insulator and 4 for graphene. The 
angular factors (/p+q.y |/ P7 ) are specific to chiral Dirac 
electron and arise in (21) as a result of summation over 
spinor components of electron wave function. 

The random phase approximation becomes exact in 
the limit of small values of dimensional parameter r s , 
defined as a ratio of characteristic Coulomb interaction 
energy to kinetic energy. For the gas of massless parti- 
cles, r s does not depend on electron density and equals 
to r s = e 2 /evF (effectively, e = (e\ + £2)/2 when two- 
dimensional electron layer is surrounded by two dielectric 
half-spaces with susceptibilities e\ and e-i). For Bi2Se3, 
r s = 0.09 with e = 40 for dielectric half-space, and ap- 
plicability of the RPA is well established (the value of r s 
for another material Bi2Te3 is close to that for Bi2Se3). 
In the case of graphene, r s is not small, but applicability 
of the RPA can be established due to smallness of the 
parameter l/g, leading to selection of bubble diagrams 
(see [16] and references therein, and also the work [26]). 

Maximal achievable amounts of doping of helical liquid 
on a surface of Bi 2 Se3 are \x ~ 0.3 eV [24], therefore at 
room temperature it can be degenerate electron liquid. 
Hence we assume that T — in all calculations below. 



IV. WAVE FUNCTION OF PLASMON 

Plasmon dispersions f2 9 , calculated numerically from 
(20)-(21) at various r s , are plotted in Fig. 1. We also 
present results for suspended graphene with r s = 8.8 (e = 
1) and for graphene embedded into SiC>2 environment 
with r s = 2.2 (e = 4). For convenience, the degeneracy 
factor g = 4 is included here into r s , since the plasmon 
dispersion in RPA (20) depend only on the combination 
gr s - 

It is seen that for small values of r s (the case of topo- 
logical insulator) the plasmon dispersion law approaches 
very closely to the upper border of the continuum of in- 
traband single-particle transitions (uj < vfq). On the 
contrary, at moderate and large r s , plasmon has well- 
defined square-root dispersion in the long-wavelength 
range. At q ~ pf, the plasmon enters the continuum 
of interband single-particle transitions (uj > 2/j — v-pq)- 
Inside the single-particle continuum, energy and mo- 
mentum conservation laws allow energy transfer between 
plasmon and single-particle excitations, so the plasmon 
acquires a finite lifetime. 

The internal structure of plasmons in Dirac electron 
gas can be investigated by looking at total weights (19) 
of intra- and interband transitions in its wave function. 




2 

a q /n 

1.5 



VI 1 1 1 1 1 1 1 1 1 1 1 1 1 

" continuum of interband excitations " 


- 
- 


y \ J? 


" / / 
J / 


* s r \ 

0.09 \^ 


1 / 




1/ y 


8.8 


L jr continuum of intraband excitations . 
t^i i i i 1 i i i i 1 i i i i 



0.5 1 2 

qiPv 

FIG. 1. Plasmon dispersions Q q at different values of r s . Con- 
tinuums of intraband (w < VFq) and interband (w > 2fi~VFq) 
single-particle excitations are shaded; pp = fi/vp is the Fermi 
momentum, n is the chemical potential. 



FIG. 3. Squared modulus of plasmon wave function |Cp q 7 | 2 
on the plane of the hole momentum p at q — 0.4pp. Top row: 
intraband channel |Cpq | at r s = 0.09 (a) and r s = 8.8 (b). 
Bottom row: interband channel |CiT| + |Cp^"| at r s = 0.09 
(c) and r s = 8.8 (d). 
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FIG. 2. Weight D++ (19) of intraband electron-hole transi- 
tions in plasmon wave function at different values of r„. 



tributions of intraband |Ci + | 2 or interband |Ci | 2 + 
|C^ + | 2 electron-hole pairs into wave function of plasmon 
with given momentum q is presented in Fig. 3. The re- 
sults are calculated for q = OApp at r s = 0.09 and 8.8. 
For other values of plasmon momentum q, the results 
are qualitatively similar. It is seen that for small values 
of r s (the case of topological insulator) the distribution 
of intraband contributions is very sharply peaked in the 
forward direction, whereas the contribution of interband 
transitions is negligible. It can be understood from the 
reason that the plasmon dispersion is very close to the 
single-particle continuum (see Fig. 1) and thus the plas- 
mon itself behave almost as single intraband electron- hole 
transition. 

At large values of parameter r s (the case of suspended 
graphene), the broad range of electron intraband tran- 
sitions in momentum space contribute to plasmon. The 
weight of interband transitions is small but not negligi- 
ble. Contributions of interband transitions form two side 
lobes, since the angular factor (/ P + q ,± |/p=p) suppresses 
interband forward scattering. 



The weight -D++ of intraband single-particle transitions 
is plotted in Fig. 2. It is seen that the undamped plas- 
mon for all values of parameter r s consists mainly of in- 
traband transitions, thus an influence of valence band on 
its properties is rather weak. When the plasmon enters 
the single-particle continuum, inter- and intraband tran- 
sitions start to contribute almost equally to its wave func- 
tion, but the plasmon undergoes strong Landau damping. 

The detailed picture of the plasmon wave function in 
momentum space C 7 7 , showing the distribution of con- 



V. CHARGE- AND SPIN-DENSITY WAVES 

When a spin-plasmon is excited in the helical liquid, 
anisotropic distribution of electron-hole pairs of the type, 
depicted in Fig. 3, arises. This distribution is shifted 
towards the plasmon momentum q. Due to the spin- 
momentum locking, the system should acquire a total 
nonzero spin polarization, perpendicular to q. A simi- 
lar situation occurs in the current-carrying state of the 
helical liquid, which turns out to be spin-polarized [9]. 
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FIG. 4. In-plane transverse spin polarization of the system (s) 
(22) in the state with single spin-plasmon excited as function 
of its momentum q at different values of r s . 



Average of spin polarization in the state |l q ) = 10) 



with a single plasmon excited is 
can be calculated using (8) as: 



(l q |er/2|l q ) and 



P77'r 
-CpqifprWlfp-y 



p+q,7' l^l/p+q^Cpq 



rn 7 



(riprj - n p+q . 7 ')- (22) 



If q is parallel to e x , only the y-component of (s) is 
nonzero. Its dependence on q for undamped plasmons 
at various r s is plotted in Fig. 4. At large enough r s , 
the spin polarization of the system is comparable to the 
whole spin of a single electron. 

Note that in the case of graphene the isospin polar- 
ization of the system appears instead of spin polariza- 
tion. Isospin polarization corresponds to nonzero aver- 
age phase shift between wave function of electrons on 
different sublatticcs. On the contrary, in pscudospin- 
unpolarizcd state this shift is zero on the average. 

In ordinary electronic system without spin-momentum 
coupling caused by spin-orbit interaction, plasmon man- 
ifests itself as charge-density wave. In helical liquid 
on surface of topological insulator, due to the spin- 
momentum locking, spin-plasmon manifests itself as cou- 
pled charge- and spin-density waves. These waves, ac- 
companying spin-plasmon with the momentum q, can 
be characterized by corresponding spatial harmonics of 
charge- (11) and spin-density 



3 q ~~ o E (•/'p+q,7'l 'l/p7 



' a p+q,7' a 



P7 



(23) 



P77' 



operators. 



To relate these quantities to the plasmon operators (8), 
we employ the unitary transformation, connecting opera- 
tors of electron-hole excitations in noninteracting system 
a p+q 7' a P7 wn; h operators of plasmons Q+. For this pur- 
pose, we also need explicit expressions for operators of 
distorted single-particle excitations in the system with 
Coulomb interaction (similarly to the work [27]). De- 
noting creation operator for such single-particle excita- 
tion with total momentum q and energy £ p + q . 7 ' — £ P7 
as J 7p q77 / j we can use the equation-of- motion method in 
the RPA and, similarly to the case of plasmons, get the 
explicit expression for it: 



U2Z7 a+ 



~ 

P'tt' 



p'+q,T' a P'l"i 



(24) 



where 



7Y77 TT £ A A 

U pp'q ~ "pp'OjtOj't' 

Vq (/p7l/p+q,7')( n P7 ~ n P+q,7') 



S 1 - V q Il(q, £p+ q , 7 ' 

(/p'+q,i"' l/p'-r) 



+q,T' 



^p+q,7' 



iS 



(25) 



The expressions (8), (16), (24) and (25) establish 
the unitary transformation at given q from opera- 
tors of electron-hole excitations in noninteracting sys- 
to operators of excitations in Coulomb- 



+ 



tern a p+qj7 ,a p7 
interacting system: plasmons and single-particle ex- 
citations ?7p q77 < ■ We can easily derive the inverse trans- 
formation of the form: 



E(^p 

P'tt' 



rr 77 

pq 



\* + 



(26) 



According to (26), we can represent the operators (11) 
and (23) of charge- and spin-density waves in the form: 



p+ = SN*U(q,n q )Q+ 

+ 
q 



Pq, 



s + = SN:n s ( qi n q )Q+ + s 



(27) 
(28) 



where the parts p"^ and 11+ are the contributions of single- 
particle excitations and are dynamically independent on 
plasmons. Here, along with the usual susceptibility (21), 
the crossed spin-density susceptibility of the helical liquid 
[11] has been introduced: 

n S (<?,w) = — ^ (/p+q,7'l/p7) 
P77' 



X (/p7l°"l/p+q,- 



'p+q.7' 



Cp7 6 



p+q,7 



iS 



(29) 



The formulas (27)-(28) show us that the average val- 
ues of and in any state with a definite number 
of plasmons vanish (similarly to the mean value of coor- 
dinate or momentum in the simplest quantum harmonic 
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oscillator). However, we can calculate their mean squares 



in the n q -plasmon state \n q 



(PqPq) = ("qlPqPql^q 



i(Q+) n ->/wr i/2 ] 

(<W+|0>, (30) 

4( S q) + l ) ( 31 ) 



««) + >^KI<«) + ^q) 



(only the in-plane transverse component s 1 - of the spin 
s is nonzero upon averaging). Since we are interested in 
plasmons only, we have subtracted the vacuum fluctua- 
tions of these quantities in the ground state |0). 

Using (27)-(28), the mean squares of charge- and spin- 
density wave amplitudes (30)-(31) can easily be calcu- 
lated: 

(Pq/£) =n cl S 2 \N ci U(q,n q )\ 2 , (32) 
(^(4) + )=n q 5 2 |iV q n^(g,0 9 )| 2 . (33) 

We can normalize the amplitudes to obtain dimension- 
less quantities A p (q) = [(/9 q p q )/n q S'p] 1 / 2 and A s (q) = 

[( s q ( s q) + )/ n q , S'p] 1 ^ 2 (p = v\ I ^ is the average electron 
density), plotted in Fig. 5. As seen, the amplitudes of 
charge- and spin-density waves are close quantitatively 
at moderate momenta and any r s . 

In [11], it was shown that, due to spin-momentum lock- 
ing, electron density and transverse component of spin 
obey an analogue of "continuity equation" . It requires 
that 

n q A p (q) = 2v F qA s (q). (34) 
Our results are in agreement with this equation. 

VI. SPIN-PLASMON SCATTERING 

Nontrivial internal structure of plasmon in Dirac elec- 
tron gas can reveal itself in a process of its scattering on 
external potentials or impurities. For the case of spin- 
plasmon, it is also interesting to consider its scattering 
on magnetic field, acting on electron spins. Since we are 
interested in spin effects here, we will consider only in- 
plane magnetic field, affecting only spins of electrons and 
not their orbital motion. 

Hamiltonian of interaction of electrons in helical liquid 
with external electric U(r) and magnetic H(r) fields are, 
respectively, 



H c = -e 2^ EV-p(/ P '7'I/p7) «p' 7 ' a P7> ( 35 ) 
pp'77' 

H m = /.jb 22 Hp/_ p (/p' 7 <|cr|/p 7 ) ap, 7 ,a p7 , (36) 
pp'77' 

where C/ q and H q are Fourier components of external 
electric and magnetic fields, fiB = e/2mc is the Bohr 
magneton. 

To calculate the probability of elastic scattering of the 
spin-plasmon with initial momentum q to the state with 
the final momentum q', we can use the Fermi golden rule, 
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FIG. 5. Normalized amplitudes of charge-density A p (a) and 
spin-density A s (b) waves in many-plasmon state as function 
of plasmons momentum q at different values of r s . 



corresponding to the Born approximation. Thus, the dif- 
ferential (with respect to the scattering angle) probabili- 
ties of this scattering on electric and magnetic fields can 
be presented in the following form: 



de 1 q q 

dw c 

~~dT 



2 l*cM)| 2 



H 



q -q 



(37) 
(38) 



where 9 is scattering angle (i.e. the angle between q and 
q'); q is absolute value of both q and q'. Here we have 
introduced electric & e (q-, &) and magnetic <f> m (q, 0) form- 
factors of spin-plasmon. 
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FIG. 6. Polar graphs of the squared modulus of electric 
|$ c (g,6>)| 2 and components \®L(q,6)\ 2 , |$i(q,6>)| 2 of mag- 
netic form-factor of spin-plasmon with the momentum q = 
0.6p P at r s = 0.09 (a) and r s = 8.8 (b). 



Using (35)-(36), we can calculate the form- factors ex- 
plicitly: 

$ e (q,e) = eJ^ r J2 «/p+q',7'l/ P +q,r) 

^ P77'r 

•r/'r 



X ^pq Cp+q'-q,q(/p+q'-q,-r|/p7/ 

x(C2jr("P7-np+q'7')» (39) 

*m(?,0) =HB\hy ((/p+q'7'kl/p+qr) 



X Cp^J ~~ Cp+q'-q,q(/p+q'-q-r| 'l/p7/ 
X (^pq')*( n P7 — n P+q'7 



). (40) 



Here Q' q = dQ q /dq is the derivative of spin-plasmon dis- 
persion law. It is convenient to project the vector 3> m 
on directions, parallel and perpendicular to the initial 
plasmon momentum q to get 4>jL and respectively. 

We consider only the form-factors, revealing the 
specifics of spin-plasmons, instead of the whole differen- 
tial probabilities (37)-(38), also dependent on the form 
of external field. In Fig. 6, angle dependencies of squared 

modulus |$cM)| 2 , |$!lM)| 2 and \$i(q,9)\ 2 of the 
form-factors at q = 0.6pf are plotted; these angular 
distributions are normalized to unity. In their calcula- 
tion, we have used the single-band approximation, since 
undamped plasmons consist mainly of interband transi- 
tions, according to Fig. 2. Results for other values of 
plasmon momentum are qualitatively the same. 

In the case of forward scattering with zero momentum 
transfer (at 9 — 0), the external electric field probes the 
total charge of plasmon, which is actually zero. Thus, 
the corresponding electric form- factor $ C (<Z, 0) = 0. Sim- 
ilarly, forward scattering on magnetic field probes the 
total spin of spin-plasmon, which is directed perpendicu- 
larly to q. Therefore, $!L(<?,0) = and $^(?,0) ^ 0. As 
for backscattering of plasmon, it is strictly prohibited, as 
for individual massless Dirac electrons [28]. 



As seen in Fig. 6, the form-factors demonstrate two 
side lobes, rather sharp at small 7' s , which can be consid- 
ered as a consequence of sharp peaking of the plasmon 
wave function (see Fig. 3). At large r s , these lobes are 
much broader. 

In the process of spin-plasmon scattering on some con- 
figurations of electric or magnetic fields, we can expect 
interplay between spatial structure of these fields and 
that of the spin-plasmon. Angular distribution of the 
scattered plasmons will incorporate both of these factors, 
according to (37)-(38). Controlling the overlap of max- 
ima of external potential and spin-plasmon form-factors, 
one can manipulate spin-plasmon scattering. 

Magnetic or nonmagnetic impurities can also create 
the external field. If the characteristic radius of the im- 
purity potential is R, the transferred momentum will be 
limited by by the order of magnitude. If qR ^ I, 
the plasmon scattering on impurities will be suppressed 
due to proximity to the regime of forward scattering. In 
the opposite limit qR <C 1, plasmons will be effectively 
scattered on considerable angles. 



VII. CONCLUSIONS 

The properties of collective excitations (plasmons or 
spin-plasmons) in two-dimensional gas of massless Dirac 
particles were studied. Two physical realizations of such 
systems were considered: electron gas in graphene and 
helical liquid on the surface of topological insulator. 
Quantum field-theoretical formalism for comprehensive 
description of spin-plasmons as composite Bose particles 
in the random phase approximation was developed. In- 
ternal structure and wave function of spin-plasmons were 
studied. 

Signatures of spin-momentum locking in helical liquid 
were considered. In particular, it was shown that exci- 
tation of a spin-plasmon induces the total nonzero spin 
polarization of the system. Moreover, coupling between 
charge- and spin-density waves, accompanying a spin- 
plasmon, was demonstrated. It was shown that ampli- 
tudes of both of these waves are close by the order of 
magnitude for spin-plasmons of intermediate momenta. 
The results of this work can be confirmed by experi- 
ments involving spin-plasmon excitation on the surface 
of topological insulator and independent measurements 
of charge and spin wave amplitudes (one of experiments 
of this type was proposed in [11]). The similar effect of 
coupling between charge and spin appears in electron gas 
with spin-orbit coupling [29], but amplitude of spin wave 
is considerable less in this case than that of charge waves 
for experimentally relevant parameters. 

Elastic scattering of spin-plasmons in helical liquid on 
electric and magnetic external fields is considered. An- 
gular distribution of scattered spin-plasmons depend on 
both the shape of the potential and the form-factor of 
the spin-plasmon, revealing its internal structure. It was 
shown that, due to the form-factor, the scattering occurs 
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into two side lobes, while forward and backward scatter- 
ing is suppressed. One can use this fact to manipulate 
spin-plasmon scattering via interplay between plasmon 
form-factor and form-factor of the external potential. It 
can be also concluded that scattering of spin-plasmons 
on long-range impurities should be very weak. 

Coupling between charge and spin waves, demon- 
strated in this article, can be used for realization of var- 
ious spintronic devices. One can perform controllable 
focusing of spin-plasmon waves and thus create regions 
with high spin polarization. Spin-polarized electrons, ac- 
cumulating in these regions, can diffuse to adjacent elec- 
trodes and be used to drive spin currents (similarly to 
[12])- 

The quantum field-theoretical approach, presented in 
this article, can be used for theoretical description of an 
influence of various external factors (impurities, external 
fields) on plasmons in Dirac electron gas. The problems 
of two-dimensional spin-plasmon optics, based on ma- 



nipulations by inhomogencitics of the system and three- 
dimensional environment, can be solved. Also the prop- 
erties of hybrid modes (plasmon polaritons [30] , plasmon- 
phonon modes [21, 22], plasmon- hole modes — plas- 
marons [31]) can be studied using this approach. 

The classical electrodynamic approach based on 
Maxwell equations and response functions cannot de- 
scribe quantum effects, arising when individual plasmons 
are emitted and adsorbed. Therefore, the quantum field- 
theoretical approach, presented in this article, should be 
especially useful for the problems of quantum plasmon- 
ics, which seems to be rather feasible in graphene-based 
structures [23]. 

The work was supported by Russian Foundation for 
Basic Research (grants 11-02-12209-ofi-m and 10-02- 
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